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Abstract Porcine deltacoronavirus (PDCoV) has been 
reported in many countries, including Hong Kong, the 
United States, South Korea, China and Thailand. In Jan¬ 
uary 2016, clinical diarrhea similar to that of porcine epi¬ 
demic diarrhea virus (PEDV) with a lower mortality rate 
was reported on a swine farm in Lao PDR. Intestine sam¬ 
ples were collected from 3-day-old pigs with clinical 
diarrhea and assayed for the presence of swine enteric 
coronaviruses. The PCR results were positive for PDCoV 
but negative for PEDV and TGEV. A phylogenetic tree 
demonstrated that PDCoV from Lao PDR was grouped 
separately from PDCoV isolates from China and the USA, 
but was more closely related to the Chinese isolates than to 
the US isolates. The full-length genome sequence of the 
novel PDCoV isolate P1_16_BTL_0116 was determined. 

Porcine deltacoronavirus (PDCoV) is an enveloped, posi¬ 
tive sense, single-stranded RNA virus belonging to the 
genus Deltacoronavirus , family Coronaviridae , order 
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Nidovirales. The family Coronaviridae is divided into four 
genera, including Alphacoroanavirus , Betacoronavirus , 
Gammacoronavirus , and Deltacoronavirus [6]. PDCoV, a 
recently described enteric pathogen, was first discovered in 
Hong Kong in 2012 [7]. In February 2014, PDCoV was 
first reported on a swine farm in Ohio, USA. Since then, the 
virus has been detected in many US states in addition to 
Asian countries, including Korea and China [2, 4, 5]. In 
this study, we determined the full-length genome sequence 
of PDCoV isolated in the Lao People’s Democratic 
Republic (Lao PDR). 

In January 2016, a clinical diarrhea outbreak was 
reported on a swine farm located in Khammouane pro¬ 
vince, Lao PDR. The herd was a farrow-to-finisher herd 
with an inventory of 1,200 sows. The breeding herd had 
four breeding-gestating rooms and six farrowing rooms. 
Farrowing houses operate an all-in/all-out system on a 
weekly basis. Clinical abnormalities including vomiting 
and watery diarrhea were observed in early January 2016 in 
a litter of 3-day-old pigs. Diarrhea was then observed in 33 
of 59 litters in the room within 12 hours. Two hundred 
fifty-seven of 603 pigs died 3 days later, and the consulting 
veterinarian ordered a planned exposure with intestines 
from pigs with diarrhea. All sows in the herd were orally 
administered infected intestine (feedback) for 3 consecu¬ 
tive days. The sow herd status resumed to normal 
approximately 3 weeks after feedback, and the breeding 
herd started to rear pigs with no clinical disease. A total of 
637 pigs were lost in the outbreak. 

Prior to feedback, 10 intestinal samples were collected 
from 3 to 5-day-old pigs displaying watery diarrhea and 
tested for viral enteric pathogens of the family Coronaviri¬ 
dae , including porcine epidemic diarrhea virus (PEDV), 
transmissible gastroenteritis virus (TGEV), and porcine 
deltacoronavirus (PDCoCV), using a PCR technique. All 


Published online: 16 July 2016 


4?) Springer 





A. Lorsirigool et al. 


intestinal samples were minced into small pieces. Total RNA 
was extracted and converted to cDNA using a Nucle- 
ospin®RNA Virus kit (Macherey-Nagel Inc., PA, USA). 
PCR was performed using specific primers for the spike 
(S) gene of PEDV [3], the nucleocapsid (N) gene of TGEV 
[1], and the membrane (M) and nucleocapsid (N) genes of 
PDCoV [5] for screening. Sequencing was performed by 
First BASE Laboratories Inc. (Selangor, Malaysia). 

The PCR results demonstrated that all samples were 
positive for PDCoV (10/10), and none were PCR positive 
for PEDV (0/10) or TGEV (0/10). After that, the full- 
length genome was sequenced using 26 overlapping 
regions of each genome (Supplementary Material). Each 
region was cloned into pGEM-T Easy Vector (Promega, 
Madison, WI, USA) and sequenced in both directions in 
triplicate as per the previously reported protocol [5]. The 5' 
and y terminal regions were determined using a kit for 
rapid amplification of 5' and 3' cDNA ends (5' and 3'- 
RACE) (Clontech, Japan). 

PDCoV designated P1_16_BTL_0116 was isolated. The 
full-length genome sequence was deposited in the Gene- 
Bank database under accession number KX118627. The 
full-length genome sequence of the isolate is 25,402 
nucleotides (nt) in length. Its genome organization resem¬ 
bles to that of all previously sequenced PDCoV genomes, 
with the gene order S'-ORFla/lb-S-envelope (E)-mem- 
brane (M)-nonstructural protein 6 (Nsp6)-nucleocapsid 


(N)-nonstructural protein 7 (Nsp7)-3 / . The lengths of 
ORFla/lb, S, E, M, and N genes are 18,786 nt, 3,477 nt, 
249 nt, 651 nt, and 1,026 nt, respectively. To determine 
the evolutionary relationships between P1_16_BTL_0116 
and other PDCoVs available in GenBank, phylogenetic 
analysis was carried out based on the full-length genome 
sequence. The results of the phylogenetic analysis 
demonstrated that the P1_16_BTL_0116 isolate belongs to 
a new group of PDCoVs that are separated from both the 
Chinese and US groups (Fig. 1). Five genes, including 
ORFla/lb, S, M, and N, were further analyzed. 

Sequence analysis demonstrated that the 
P1_16_BTL_0116 isolate is closely related to isolates from 
China, with 97.1-97.7 % and 93.7-95.0 % sequence iden¬ 
tity at the nucleotide and amino acid level, respectively 
(Table 1). Moreover, the P1_16_BTL_0116 isolate is clo¬ 
sely related to isolates from the USA, with 97.2-97.3 % 
and 94.0-94.1 % sequence identity at the nucleotide and 
amino acid level, respectively. 

Based on the ORFla/lb gene, pairwise nucleotide and 
amino acid sequence identity values between the 
P1_16_BTL_0116 isolate and isolates in the China-like 
group were in the range of 97.2-97.9 % and 95.8-96.7 %, 
respectively, and and those with the US-like group were 
913-91A % and 96.0-96.1 %, respectively. Moreover, 84 
substitutions and five deletions of two ( 401 LK 402 ) and three 
( 758 PVG 760 ) amino acids at positions 401 to 402 and 758 to 
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Fig. 1 Phylogenetic analysis based on the nucleotide sequence of the 
full-length genome of porcine deltacoronavirus (PDCoV), performed 
using the Bayesian Markov chain Monte Carlo (BMCMC) method. 
The P1_16_BTL_0116 isolate, the isolate from Lao PDR, is 


highlighted in yellow. The red and blue underlining indicates the 
US and China PDCoV group, respectively. The reference sequences 
obtained from GenBank are indicated by strain name and accession 
number (color figure online) 
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Genome sequence of porcine deltacoronavirus from Lao PDR 


Table 1 Comparison of the nucleotide and amino acid sequence identities (%) based on the full-length genome sequence of the 
P1_16_BTL_0116 isolate with those of porcine deltacoronavirus (PDCoV) isolates previously reported from the USA and China 


PDCoV isolate 

Gene 

US PDCoV 


China PDCoV 


Nucleotide (%) 

Amino acid (%) 

Nucleotide (%) 

Amino acid (%) 

P 1_ 16_BTL_0116 

Full-length genome 

97.2-97.3 

94.0-94.1 

97.1-97.7 

93.7-95.0 


ORFla/lb 

97.3-97.4 

96.0-96.1 

97.2-97.9 

95.8-96.7 


S 

96.2-96.5 

98.2-99.1 

95.7-96.9 

98.5-99.1 


E 

85.6-98.4 

85.5-97.6 

85.2-98.0 

85.3-97.6 


M 

98.0-98.3 

99.5 

98.1-98.7 

99.5 


N 

97.8-98.1 

98.2-99.1 

98.1-98.7 

98.5-99.1 


760 were observed between P1_16_BTL_0116 and isolates 
in the China-like and US-like groups. 

Based on the S gene, the P1_16_BTL_0116 isolate is clo¬ 
sely related to Chinese PDCoVs, with 95.7-96.9 % and 98.5- 
99.1 % nucleotide and amino acid sequence identity, 
respectively. It also shared sequence similarity with the US 
PDCoV isolate, with 96.2-96.5 % and 98.2-99.1 % identity at 
the nucleotide and amino acid level, respectively. Twenty - 
three and 25 amino acid substitutions were observed in the 
Laos PDCoV isolate when compared to Chinese and US 
PDCoV isolates, respectively. In addition, the 
P1_16_BTL_0116 isolate contained a deletion of one amino 
acid (N) at position 51, similar to isolates in the China-like 
group. 

In the M gene, compared to Chinese PDCoV isolates, the 
Lao PDCoV isolate had 98.1-98.7 % and 99.5 % nucleotide 
and amino acid sequence identity, respectively. The Laos 
PDCoV isolate had 98.0-98.3 % nucleotide and 99.5 amino 
acid sequence identity to the US PDCoV. Only one amino 
acid substitution was observed between the Laos PDCoV 
isolate and the Chinese and US PDCoV isolates. 

In the N gene, the nucleotide and amino acid sequences of 
the Lao PDCoV isolate were 98.1-98.7 % and 98.5-99.1 % 
identical, respectively, to those of Chinese PDCoV isolates 
and 97.8-98.1 % and 98.2-99.1 % identical, respectively, to 
those of the US PDCoV isolates. Four amino acid substitu¬ 
tions were observed between the Laos PDCoV isolate and the 
Chinese and US PDCoV isolates. 

In this study, the full-length genome sequence of a PDCoV 
isolate from Lao DPR was determined. The findings suggest 
the presence of PDCoV in Lao PDR causing a diarrhea out¬ 
break similar to that of PEDV. Phylogentic analysis showed 
that this novel isolate cluster separately from Chinese and US 
PDCoV isolates but was more closely related to the Chinese 
PDCoV that to the US isolates. Therefore, there is an urgent 
need to monitor the infection status of PDCoV and the 
molecular epidemiology in Lao PDR. 

Nucleotide sequence accession number The complete 
genome sequences of P1_16_BTL_0116 isolate has been 
deposited in GenBank under accession number KX118627. 
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